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ABSTRACT 

We  continue  our  investigation  of  overcoming  Gibbs  phenomenon,  i.e.,  to  obtain  exponen¬ 
tial  accuracy  at  all  points  (including  at  the  discontinuities  themselves),  from  the  knowledge 
of  a  spectral  partial  sum  of  a  discontinuous  but  piecewise  analytic  function.  We  show  that  if 
we  are  given  the  first  N  Gegenbauer  expansion  coefficients,  based  on  the  Gegenbauer  polyno¬ 
mials  Cj*(x)  with  the  weight  function  (1  —  for  any  constant  p  >  0,  of  an  L\  function 

/(x),  we  can  construct  an  exponentially  convergent  approximation  to  the  point  values  of 
f{x)  in  any  sub-interval  in  which  the  function  is  analytic.  The  proof  covers  the  cases  of 
Chebyshev  or  Legendre  partial  sums,  which  are  most  common  in  applications. 
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1  Introduction 


In  this  paper  we  continue  our  investigation  of  overcoming  the  Gibbs  phenomenon,  i.e.,  recov¬ 
ering  pointwise  exponential  accuracy  at  all  points  including  at  the  discontinuities  themselves, 
from  the  knowledge  of  a  spectral  partial  sum  of-  a  discontinuous  but  piecewise  analytic  func¬ 
tion,  which  we  started  in  [4],  [5],  and  [6]. 

Spectral  approximations,  such  as  the  Fourier  approximation  based  upon  trigonometric 
polynomials  for  periodic  problems,  and  the  Chebyshev,  Legendre  or  the  general  Gegenbauer 
approximation  based  upon  polynomials  for  non-periodic  problems,  are  exponentially  accurate 
for  analytic  functions  [3],  [2].  However,  for  discontinuous  but  piecewise  analytic  functions, 
the  spectral  partial  sum  approximates  the  function  poorly  throughout  the  domain.  Away 
from  the  discontinuity  only  first  order  accuracy  is  achieved.  Near  the  discontinuity  there  are 
0(1)  oscillations  which  do  not  decrease  with  N,  the  number  of  terms  retained  in  the  spectral 
sum.  This  is  known  as  the  Gibbs  phenomenon. 

Our  framework  in  [4],  [5]  and  [6]  to  overcome  Gibbs  phenomenon  to  obtain  exponential 
accuracy  at  all  points  for  piecewise  analytic  functions  relies  heavily  on  using  the  Gegenbauer 
polynomials  0^(x),  which  are  orthogonal  in  [—1, 1]  with  the  weight  function  (1  —  ,  for 

large  A.  We  assume  that  the  first  —N  <  k  <  N  Fourier  coefficients,  or  the  first  0  <  k  <  N 
Legendre  coefficients,  of  a  discontinuous  but  piecewise  analytic  function,  are  given.  The 
procedure  consists  of  two  steps: 


1.  Using  the  given  spectral  partial  sum  of  to  recover  the  first  m  ~  Gegenbauer 
expansion  coefficients,  based  on  a  sub-interval  [a,  6]  C  [— U  1]  in  which  the  function 
is  presumably  analytic,  with  exponential  accuracy.  This  can  be  achieved  for 
any  Li  function,  as  long  as  we  choose  A  in  the  weight  function  of  Gegenbauer 
polynomials  to  be  proportional  to  N.  The  error  incurred  at  this  stage  is  called 
the  truncation  error, 

2.  For  an  analytic  function  in  [a,  6],  proving  the  exponential  convergence  of  its  Gegen- 

JtBt  j 

'  I 


□ 

□ 


-ides 

or 


bauer  expansion,  when  the  parameter  A  in  the  weight  function  is  proportional  to 
the  number  of  terms  retained  in  the  expansion.  The  error  at  this  stage  is  labeled 
the  regularization  error. 

In  [6]  we  demonstrated  this  procedure  in  the  case  of  a  discontinuous  but  piecewise  analytic 
function,  provided  its  Fourier  or  Legendre  spectral  partial  sum  is  given. 

The  proof  of  the  Legendre  case  in  [6]  is  based  upon  first  expanding  the  Legendre  poly¬ 
nomial  Pk{x)  =  Ck{x)  into  its  Fourier  series.  It  was  essential  in  this  proof  that  the  Fourier 
expansion  for  the  Legendre  polynomial  Ps{x),  for  large  N,  contains  lower  terms  that  decay 
exponentially  with  N  (formula  (2.13)  in  [6]).  Unfortunately,  it  seems  that  this  fact  is  true 
only  for  Legendre  polynomials,  probably  because  their  weight  function  is  special  (=  1).  It 
seems  not  true  for  other  Gegenbauer  polynomials,  such  as  Chebyshev  polynomials.  In  an 
earlier  version  of  [6],  we  quoted  a  formula  (7.354,  page  836  of  [7])  to  this  effect  for  Chebyshev 
polynomials.  However,  it  is  doubtful  that  Formula  7.354  of  [7]  is  correct. 

In  this  paper,  we  will  consider  the  case  of  general  Gegenbauer  spectral  methods,  with 
Chebyshev  and  Legendre  methods  as  special  cases.  We  assume  that  f{x)  is  an  L\  function 
on  [—1, 1]  and  analytic  in  a  subinterval  [a,  6]  C  [— 1,  !]•  We  also  assume  that  the  Gegenbauer 
partial  sum  of  /(x),  based  upon  the  Gegenbauer  polynomials  C^{x)  with  the  weight  function 
(1  —  for  any  constant  ft  >  0,  over  the  full  interval  [—1, 1],  is  known.  The  objective  is 

to  recover  exponentially  accurate  point  values  over  the  subinterval  [a,  6]  of  analyticity. 

We  will  follow  the  same  path  as  in  [6].  Basically  we  will  show  that  the  first  0  <  k  <  N 
Gegenbauer  expansion  coefficients,  based  on  the  Gegenbauer  polynomials  C^(x)  for  any 
constant  /i  >  0,  contain  enough  information,  such  that  a  different,  rapidly  converging 
Gegenbauer  expansion  in  the  subinterval  [o,6],  with  the  parameter  A  in  the  weight  func¬ 
tion  (1  —  being  proportional  to  N,  can  be  constructed.  As  before,  we  will  separate 

the  analysis  of  the  error  into  two  parts:  truncation  error  and  regularization  error.  Trunca¬ 
tion  error  measures  the  difference  between  the  exact  Gegenbauer  coefficients  with  A  ~  yV, 
and  those  obtained  by  using  the  spectral  partial  sum.  This  will  be  investigated  in  Section  3. 
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The  regularization  error  measures  the  difference  between  the  Gegenbauer  expansion  using 
the  first  few  Gegenbauer  coefficients  with  X  N,  and  the  function  itself  in  a  sub-interval 
[a,  6],  in  which  the  function  is  assumed  analytic.  This  error  is  estimated  in  [6]  and  we  will 
simply  quote  the  result  in  Section  4.  The  results  are  summarized  in  Theorem  4.3,  and  some 
remarks  are  also  given  in  Section  4.  Section  5  contains  two  numerical  examples  to  illustrate 
our  results.  In  Section  2  we  collect  some  useful  properties  of  Gegenbauer  polynomials  to  be 
used  later. 

Throughout  this  paper,  we  will  use  A  to  denote  a  generic  constant  or  at  most  a  polynomial 
in  the  growing  parameters,  as  will  be  indicated  in  the  text.  It  may  not  be  the  same  at  different 
locations. 

2  Preliminaries 


In  this  section  we  collect  some  useful  results  about  the  Gegenbauer  polynomials,  to  be  used 
in  later  sections.  We  rely  heavily  on  the  standardization  in  Bateman  [1]. 

Definition  2.1.  The  Gegenbauer  polynomial  is  defined  by 

(1  -  =  G(A,n)^  [(1  -  (2.1) 

where  G(A,  n)  is  given  by 


G(A,n) 


(-l)’*r(A-|-i)r(n-|-2A) 
2’‘n!r(2A)r(n  +  A  ^ 


(2.2) 

□ 


Formula  (2.1)  is  also  called  the  Rodrigues’  formula  [2,  page  175]. 
Under  this  definition  we  have 


r(n  -1-  2A) 
n!r(2A) 


and 

IG.V)!  <  Gi(l),  -1<X<1 


(2.3) 

(2.4) 
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The  Gegenbauer  polynomials  are  orthogonal  under  their  weight  function  (1  — 


jjl  -  =  S,,„kt 


where 


ht  = 


r(A  +  i) 


-nV‘,r{A)(n  +  A) 

We  will  need  to  use  heavily  the  asymptotics  of  the  Gegenbauer  polynomials  for  large  n 
and  A.  For  this  we  need  the  well-known  Stirling’s  formula: 


(27r)^a:*+ie"*  <  r(x  +  1)  <  (27r)* x  >  1 


mma  2.2.  There  exists  a  constant  A  independent  of  A  and  n  such  that 

The  proof  follows  from  (2.6)  and  the  Stirling’s  formula  (2.7). 


We  also  need  the  following  lemma,  which  is  easily  obtained  from  the  Rodrigues’  formula 


(2.1): 


8.3.  For  any  A  >  1  we  have: 


^[(l-xY-^C„^(x)]  = 


G’(A,n) 

G(A- 1,71-1-1) 


(i-xY”C„^;’(x) 


The  proof  follows  from  taking  one  derivative  ^  on  both  sides  of  the  Rodrigues’  formula 
(2.1),  and  then  using  it  again  on  the  right  hand  side. 


Finally,  we  would  need  to  use  the  following  formula  [2,  page  176]: 


(2.10) 
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3  Truncation  Error  in  a  Sub-interval 


Consider  an  arbitrary  L\  function  /(x)  defined  in  [—1, 1].  Suppose  that  the  first  0  <  k  <  N 
Gegenbauer  coefficients,  based  upon  the  Gegenbauer  polynomials  Cl^(x)  with  the  weight 
function  (1  —  for  any  constant  p  >  0,  over  the  full  interval  [—1, 1],  are  given: 

r(fc)  =  ~  Jjl  -  x^Y~hc^{x)f{x)dx,  0  <  fc  <  iV  (3.1) 

We  are  interested  in  finding  the  Gegenbauer  expansion  of  /(x),  with  A  ~  A^,  based  on  a 
sub-interval  [a,  6]  C  [— 1, 1].  We  start  by  introducing  the  local  variable 
Definition  3.1.  The  local  variable  ^  is  defined  by 

X  =  x(0  =  ^ 

where 

b  —  a  .  b  +  a 
<  =  —  «  =  — 

Thus  when  n<x<6,  — 


(3.2) 

(3.3) 

□ 


We  consider  functions  /(x)  satisfying 

A 

Assumption  3.2.  |/'*(*)|  S  A  independent  of  k. 

□ 

We  remark  that  if  /(x)  is  an  Li  function  this  assumption  is  fulfilled. 

Since  we  know  the  first  yV  +  1  Gegenbauer  coefficients,  f'*{k)  for  0  <  k  <  we  define 
the  Gegenbauer  partial  sum: 

/Wx)  =  i;/”(*:)cnx)  (3.4) 

Note  that  /jv(x)  does  not  converge  fast  to  /(x)  if  there  exist  discontinuities  inside  the  domain. 


The  function  /(x)  has  also  a  Gegenbauer  expansion  in  a  sub-interval  [a,  6],  with  A  ~  ^V. 
With  f  and  8  defined  in  (3.2)-(3.3),  we  have 

+  ^)  =  -1  <  ^  <  1  (3.5) 

/sO 

A  - 

where  the  Gegenbauer  coefficients  /,  (/)  are  defined  by 

t(t) = ^  fji  + m  (s.e) 

Of  course,  we  do  not  have  f^(l)  at  our  disposal,  but  only  an  approximation  based  on  the 
Gegenbauer  parial  sum  thus  we  have 

^  X'/i  -  + m  (:!.7) 

How  well  do  g^[l)  approximate  To  answer  this  question  we  define 

Definition  3.3.  The  truncation  error  is  defined  by 

m 

TE(Km,N,,)=  -  m)Ct{0\  (3-8) 

"  i=0 

where  f^{l)  are  defined  by  (3.6)  and  g^{l)  are  defined  by  (3.7). 

□ 


The  truncation  error  is  the  measure  of  the  distance  between  the  true  Gegenbauer  ex¬ 
pansion  in  the  interval  [a,  6]  and  its  approximation  based  on  the  Gegenbauer  partial  sum  in 
1-1,1]. 

We  first  have  the  following  lemma: 

Lemma  3.4.  The  truncation  error  can  be  estimated  by 

|Ci'(i) 


TE(X,m,N,t)<  x:  irwii; 

g=yv+i  1=0 

Proof;  FVom  (3.6)  and  (3.7)  we  have 

1 


(3.9) 


f,Hi)  -  9X‘)  =  ^  /,(!  -  + <)  -  «(<« + sm 


(3.10) 


Substituting  (3.10)  into  (3.8),  recalling  (2.4)  and 


m  +  +  nq)C^(t^  +  S) 

9=N+1 


we  obtain  (3.9). 


For  simplicity  of  notations  we  denote: 

F,"''  =  X',(l  -  (e{  +  S)d( 

In  order  to  estimate  this  term  we  start  with  the  following 
Lemma  3.5.  If  we  denote 

pKi 


G(A,0 

where  G(A,  /)  is  defined  by  (2.2),  then  we  have  the  following  recursive  formula: 

1 


jx,i  ^ 

’  2(<7  +  /^)€ 


[,A_-U+1  _  _  A  >  1,  ,  >  1 


(3.11) 


□ 


(3.12) 


(3.13) 


(3.14) 


Proof:  By  the  dehoitioo  of  in  (3.13)-(3.12),  we  have 

=  ■  2ut^  /->  -  + *)] 

=  I !('  -  «>1  + *)1  ^ 

=  2W-1./Vi)(,T^  [C^.(ef  +  «)  -  +  «)]  ^ 

_  ^  r  rA-l./+l  rA-l,/+ll 

-  2(9  +  ^)cI^’-’  J 

where  we  have  used  (2.10)  for  the  second  equality;  integration  by  parts  for  the  third  equality 
(the  boundary  terms  vanish  because  of  the  term  (1  —  ^*)‘'“»  with  A  >  1);  formula  (2.9)  for 
the  fourth  equality;  and  the  definition  (3.12)-(3.13)  for  the  last  equality. 


□ 
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We  can  now  obtain  the  following  recursive  estimate  for 
Lemma  3.6.  The  7^’*  defined  by  (3.13)  satisfies  the  following  estimate 

|/M  <  ^  ^  max  L  i<miE 

I  «  I  fi 9-j<p<9+i  I  ^  ' 

Proof:  We  use  induction  on  j.  The  estimate  is  clearly  valid  for  j  —  ^ 


kM  ^  ir(  »  j  ^  min(A,9)  (3.15) 

I  »  I  f^}Y[^q n 9-J<P<9+J  '  ^  ' 

Proof:  We  use  induction  on  j.  The  estimate  is  clearly  valid  for  j  —  0.  Assume  that  it  is 
valid  for  j  =  jo  <  min(A,g)  —  1,  then 

|/.,,|  <  r(,+y-n-io)  |,A-ij+»| 

1  *  I  C'»r(^  +  /i+l)  9-io<P<9+JO  •  ^  I 

^  ^(9  4*  ~b  1  jo)  ninv  ^ _  [ /A-jo-li*+io+l  _  r-^-Jo-l.^+io+l] 

~  C'»r(9  +  /i+l)  9-io<P<9+io  2(p  + /i)c  ^  •* 

^  ^(9  "b  p  ~b  1  Jo)  _ 1 _  fl /■^— io-l.^+io+i  I  I  I /■^-Jo-i.i+jo+l  I 

~  c^r(^  +  p  +  l)  2(9  —  ^0  +  p)e  9-io<p<9+io  I 

"b  /*  ~  Jo)  mav  Ir-'-is-l.^+JO+ll 


<  "b  /*  Jo)  max  l/-'-i9-J>^+>o+i| 

—  eA)+ir(9  + /i  +  1)  9-jo-i<p<9+io+i  •  I 

where  we  have  used  (3.14)  for  the  second  inequality.  All  other  steps  are  simple  inequalities. 
This  finishes  the  induction. 


Frcnn  the  previous  lenuna  we  can  get  the  following  estimate: 
mma  3.7.  For  defined  in  (3.12)  with  A  <  q  we  have  the  following  estimate; 


lpA./|^/(g-A)  |g(A,OI 

1  ’  1  -  c-^^q)  ■  |(7(0,/  +  A 


(3.16) 


where  A  grows  at  most  as  q***"^. 

Proof:  For  simplicity  and  without  loss  of  generality  we  assume  A  is  an  integer.  Since  A  <  q, 
we  can  take  j  =  A  in  (3.15)  to  arrive  at 


UgA  <  ^  ^  ^ max  < 

\  1  I  9-A<p<9+A  I  ^  I 


I  »  I  9-A<P<9+'^  '  ^  '  C'^r'(q)  9- 

By  the  definition  (3.12)-(3.13),  we  have,  for  q  —  A  <  p  <  q  +  A, 


C'^r(q)  9--^<P<9+A  I  ^  I 
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-  \GiO,U  \)\  J-i'^  ^ 


<  liLLy.Hp+M — 1 _ Pii-er^d^ 

~  (/  +  A)!  p\r{2(jL)  \G{0,U\)\J.x'^  ^ 

<  /I— L_ 

~  |G(0,/  +  A)| 

where  for  the  second  inequality  we  have  used  (2.4)  and  for  the  third  inequality  we  have  used 
(2.3).  Clearly  /4  is  a  constant  if  |  and  A  grows  at  most  as  \{  p  >  \.  Invoking 
(3.13)  again  we  obtain  (3.16). 


Using  Stirling’s  formula  we  can  now  easily  get; 

Lemma  3.8.  For  I  <in<  N  and  q  >  N,vie  have 

If  AJ|  <  (m  +  2Ar-^>  l_ 

I  ’  I  -  (2fA)Am” 

where  A  again  grows  at  most  as  (m  +  A)^^^'*"’. 

Proof;  Starting  from  (3.16)  and  using  the  definition  (2.2),  we  obtain; 


(3.17) 


\px,i\  <  .r(<y-A)  |G(A,/)| 

I  »  I  -  e^r(q)  \G{0,l  +  X 


<  r(A  +  |)r(f  +  2A)  2'+^(/  +  A)!r(/  +  A-f  I) 

c>r(q)  2W(2A)r(/  +  A  +  i)  r(/  +  A) 


<  A 


r{q-x)  r(A)r(/ +  2A)2'' 


c^r(9)  Z!r(2A) 

^  r(g-A)  r(A)r(m  +  2A)2^ 

-  c^r(9)  ■  m!r(2A) 

<  A-^e-^(m  +  2A)’"+^^e-(’"+^^)2^ 

~  m’"e“”*(2A)*^e“^^ 

(m  +  2A)’"-^^^  _L 

~  (2cA)^m"*  q^ 

where  we  have  used  (2.2)  in  the  second  inequality;  the  monotonicity  with  respect  to  I  in  the 
fourth  ineqmdity;  and  the  Stirling’s  formula  (2.7)  for  the  fifth  inequality. 
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We  are  now  ready  for  the  main  theorem  of  this  section; 

Theorem  3.9.  Let  the  truncation  error  be  defined  in  (3.8).  Let  X  =  aeN  and  m  =  ^(N 
with  0  <  a, /?  <  1,  then 

(3.18) 

where  A  grows  at  most  as  In  particular,  if  a  =  ^  <  ^,  then 

TE{a€N,a€N,N,e)<Aq^^  (3.19) 


where 


Proof;  The  theorem  follows  from  (3.9),  the  Assumption  3.2,  (2.8),  and  (3.17). 


(3.20) 


□ 


4  Regularization  Error  and  the  Main  Theorem 


The  second  part  of  the  error,  which  is  called  tfie  regularization  error  and  is  caused  by  using 
a  finite  Gegenbauer  expansion  based  on  a  sub-interval  [a,  6]  C  [—1,1],  to  approximate  a 
function  f{x)  which  is  assumed  analytic  in  this  sub-interval,  has  been  studied  in  [6].  We  will 
thus  just  qhote  the  result. 

We  assume  that  f{x)  is  an  analytic  function  on  [a,  6]  satisfying 
Assumption  4.1.  There  exists  constants  p  >  I  and  C{p)  such  that,  for  every  A:  >  0, 

(4.1) 
P 

□ 


max 

o<®<6 


This  is  a  standard  assumption  for  analytic  functions,  p  is  the  distance  from  [a,  6]  to  the 
nearest  singularity  of  f{x)  in  the  complex  plane  (see  for  example  [8]).  Let  us  consider  the 
Gegenbauer  partial  sum  of  the  first  m  terms  for  the  function  /(c^  -f-  ^): 

m 

1=0 
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with  e  and  S  defined  by  (3.2)  and  (3.3),  and  the  Gegenbauer  coefficients  based  on  [a,b] 
defined  by 

f.Hl)  =  ^ /_' (1  -«’)*-'/(<«  +  S)Cl'U)di  (4.3) 

The  regularization  error  in  the  maximum  norm  is  defined  by; 


REiX,  m,  e)  =  m^  \f{e(  +  6)  -  /,^(/)C,^(0 

■  I  1=0 


(4.4) 


We  have  the  following  result  for  the  estimation  of  the  regularization  error,  when  A  ~  rn 

[6]: 

Theorem  4.2.  Assume  A  =  7m  where  7  is  a  positive  constant.  If  f{x)  is  analytic  in 
[0,6]  C  [—1, 1]  satisfying  the  Assumption  4.1,  then  the  regularization  error  defined  in  (4.4) 
can  be  bounded  by 

RE(‘ym^m,e)  <  Aq”'  (4.5) 


where  q  is  given  by 


which  is  always  less  than  1. 
constant,  then 


1  +7)'+^  '  ' 

In  particular,  if  7  =  1  and  m  =  0eN  where  ^  is  a  positive 


RE(ffN,ffN,c)  <  Ai,'" 


(4.7) 


with 


(4.8) 

D 


We  can  now  combine  the  estimates  for  truncation  errors  and  regularization  errors  to 
obtain  the  following  main  theorem  of  this  paper: 

Theorem  4.3.  (Removal  of  the  Gibbs  Phenomenon  for  the  sub-interval  case  of  Gegenbauer 
partial  sum). 
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Consider  a  Li  function  f{x)  on  [—1,1],  which  is  analytic  in  a  sub-interval  [a,  6]  C  [—1,1] 
and  satisfies  Assumption  4.1.  Assume  that  the  first  N  +  1  Gegenbauer  coefficients 

L 

for  p  >  0,  axe  known.  Let  0  <  /  <  m  be  the  Gegenbauer  expansion  coefficients,  defined 
in  (3.7),  based  on  the  sub-interval  [a,  6],  of  the  Gegenbauer  partial  sum  /^(x)  in  (3.4).  Then 


for  A  =  m  =  lieN  with  <  ^,  we  have 


max 


1=0 


<  ^  (4"  +  Ir) 


(4.9) 


where 


<1, 


-(F' 


and  A  grows  at  most  as 

Proof:  Just  combine  the  results  of  Theorems  3.9  and  4.2. 


‘i2pj 


<  1 


□ 


We  now  give  two  remarks: 

Remark  4.3.  Comparing  with  the  Legendre  case  in  [6],  we  can  see  that  the  current  proof  is 
less  sharp  (missing  a  factor  of  A  in  the  truncation  error  qr)-  The  main  loss  in  this  sharpness 
is  in  the  estimate  (3.15). 


Remark  4.4.  No  attempt  has  been  made  to  optimize  the  parameters. 


5  Numerical  Results 


□ 


□ 


In  this  section  we  give  two  numerical  examples  to  illustrate  our  result.  We  will  test  Chebyshev 
series  because  these  are  used  most  often  in  practice.  Notice  that  the  Chebyshev  polynomials 
are  just  Gegenbauer  polynomials  with  p  =  0  module  a  constant:  Tk{x)  =  |C'°(x). 
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<  6 


Example  S.l.  We  take  the  simple  step  function 

^  '  [0,  otherwise 

and  assume  that  we  know  the  first  iV  +  1  Chebyshev  coefficients  of  f{x): 

P>{k)  =  —  [\\-  x^)-'^Tk{x)f{x)dx,  0<k<N 

TTCt  J-1 


(5.1) 


(5.2) 


where 


Ck 


-M 


if 

if 


We  then  form  the  Chebyshev  partial  sum 


^•  =  0 
k  >  1 


(5.3) 


rtw=f:/"wn(*)  (5.4) 

k=0 

and  then  compute  the  approximate  Gegenbauer  expansion  coefficient  based  on  the  sub¬ 
interval  [o,  6]  defined  by  (3.7): 


g^W  =  ^  /.' 0  -  +  s)d(  (5.5) 

With  these  Gegenbauer  coefficients,  we  can  finally  compute  the  uniformly  accurate  ap¬ 
proximation  on  [a,  6]  defined  by 

m 

gM  =  E  (5-6) 

t=0 

Numerical  experiments  (for  various  functions)  seem  to  indicate  that 

m  =  0.1fN,  A  =  0.2eN  (5.7) 

are  good  choices.  Notice  that  in  our  proof  we  did  not  attempt  to  optimize  these  parameters. 
For  consistency  we  will  use  (5.7)  for  both  examples. 

For  this  special  function  (5.1),  there  is  no  regularization  error.  Hence  all  we  see  is 
the  truncation  error.  In  Fig.  1,  left,  we  show  the  errors  of  a  middle  sub-interval  [a,  ft]  = 
[—0.5, 0.5],  and  in  Fig.  1,  right,  we  show  that  of  a  one-sided  sub-interval  [a,  ft]  =  [0, 1].  We 
can  clearly  see  good  convergence  for  both  cases. 


13 


Fig.  1;  Errors  in  log  scale,  f{x)  defind  by  [a,  i]  =  [—0.5, 0.5]  (left)  and  [a,  5]  =  [0, 1] 

(right).  A  =  0.2clV  and  m  =  0.1  clV.  N  =  20,40,80, 160. 


Since  there  is  no  regularization  error  for  this  example,  and  the  truncation  error  is  smaller 
for  small  m,  we  also  plot  the  errors  for  m  =  1  and  A  =  0.2£iV  in  Fig.  2.  We  can  see 
that  the  errors  are  now  much  smaller  than  those  in  Fig.  1.  C  '  'ourse  for  general  functions 
regularization  errors  must  balance  with  truncation  errors,  so  we  cannot  expect  m  =  1  to 


work  for  the  general  case. 


Fig.  2:  Errors  in  log  scale,  /(x)  defind  by  (5.1).  [a,  6]  =  [—0.5, 0.5]  (left)  and  [a,  5]  =  [0, 1] 
(right).  A  =  0.2cAr  and  m  =  1.  N  =  20,40,80, 160. 

Example  5,2.  In  the  second  example  we  take  the  the  following  function 


/( 


■'•(f 


.i(cos(x)),  if  a  <  X  <  b 
otherwise 


(5.8) 


Again  we  assume  that  we  know  the  first  N  +  l  Chebyshev  coefficients  of  /(z)  defined  by 
(5.2). 


This  time  both  truncation  error  and  regularization  error  exist.  We  again  pick  two  cases 
with  middle  as  well  as  one-sided  sub-intervals.  FVom  Fig.  3  we  can  see  similar  results  as  in 
the  previous  example,  Fig.  1 . 


Fig.  3:  Errors  in  log  scale,  /(z)  defind  by  (5.8).  [a,  6]  =  [—0.5, 0.5]  (left)  and  [a,  6]  =  [0,1] 
(right).  A  =  0.2eyV  and  m  =  1.  N  =  20,40,80,160. 

These  examples  illustrate  well  the  good  convergence  behavior  of  our  approach. 

6  Concluding  Remarks 

We  have  proven  the  exponential  convergence  in  the  maximum  norm,  of  a  reconstruction 
procedure  using  Gegenbauer  series  based  on  C^{x)  with  large  A,  for  any  Li  function  in  any 
sub-interval  [a,  6]  in  which  the  function  is  analytic,  if  we  are  given  the  first  N  expansion 
coefficients  of  this  function  over  the  full  interval  [—1,1]  based  on  Chebyshev,  Legendre,  or 
any  other  Gegenbauer  polynomial  basis.  Numerical  examples  are  also  given. 
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